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Abstract 

We investigate the structure of dynamics of large self-gravitating 
astrophysical systems using a self-interacting two-component model. 
We consider two cases, galaxy clusters and cosmic walls, for illustra- 
tions. In both cases stability analyses are conducted using perturbative 
expansion. We have found that waves and solitons are easily gener- 
ated in these systems. Our analysis shows that dark matter can be 
Jeans unstable in the very inner regions of galaxy clusters if it has a 
large internal degree of freedom. The dark matter core may collapse 
under external perturbations. We also discuss dark-matter oscillations 
in galaxy clusters and how mode growth and decay lead to heating 
of intracluster medium. Our analysis shows that dark-matter solitons 
with both positive and negative amplitudes can be excited in cosmic 
walls. Resonances in soliton interaction could enhance gas condensa- 
tion. The co-existence of the two types of dark-matter solitons implies 
that bright filaments can arise in dark voids. 
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1 Multi-component Self-gravitating Systems 



The Universe is a hierarchy of self-gravitating multi-component systems: 
from walls, super clusters, galaxy clusters and groups to galaxies, star clus- 
ters, stars and sub-stellar objects. The dark energy provides a repulsive pres- 
sure on very large scales, but dark matter and baryons adhere gravitationally 
to form smaller structures. Single-component systems and multi-component 
systems possess very different properties because of the complexities in dy- 
namical interactions between the components, and between the components 
and the environment. For instance, dark matter does not radiate electro- 
magnetic waves but participates in gravitational interaction with baryons 
(gas). However, the baryonic gas, which has a much smaller inertia than the 
dark matter, regulates the energy exchange of the system with the environ- 
ment. We are still unsure as to the nature of the dark matter and how many 
dark-matter species are present. We are also not certain if dark matter is 
strongly self-interacting or essentially collisionless. 

Here, we report the recent finding of our study of two example multi- 
component self-gravitating systems: galaxy clusters and cosmic walls. We 
show the observational characteristics of galaxy clusters in our two-compoent 
model with self-interacting dark-matter. We also demonstrate the formation 
of solitary waves in cosmic walls when dark matter is self-interacting. 



2 Galaxy Clusters 

2.1 Self-interacting Multi-component Model 

Galaxy clusters consist of about 85% dark matter and 10% hot ionised gas 
in mass (see e.g. Vikhlinin et al. 2006). The trapped baryons in stars and 
galaxies have an insignificant mass contribution, and thus play a less impor- 
tant role in the cluster dynamics. For simplicity. If we ignore the magnetic 
fields, cosmic rays and galactic outflow, galaxy clusters can be treated as 
two-component systems consisting only of dark matter and baryons. The 
cluster structure and dynamics is then determined by the conservation equa- 
tions: 
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(i = 1, 2 for gas and dark matter respectively), where Vi the bulk velocity and 

Ci the energy gain/loss rate. For self-interacting dark matter and gas, the 
effective pressure pi, density pi, velocity dispersion af and specific entropy 
Si are related by pi = piaf = Sip^' , where 7i is the adiabatic index and the 
degree of freedom Fi = 2/(7i — 1). The gravitational force f = — V^, where 
the potential ^ satisfies the Poisson Equation. 

For stationary clusters, d/dt — >■ 0, and the radial structure profiles of 
the cluster is determined by 
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where G is the gravitational constant (see Saxton and Wu 2008). No ra- 
diative loss from dark matter implies C2 = For the radiative loss from 
gas, a parametrisation Ci = Apf[a\)^ is appropriate. Assume that free-free 
emission is the dominant cooling process. Then we have a = 2 and (3 = 1/2, 
and A is set by the chemical abundances. The presence of radiative cooling 
of gas causes a pressure depletion in the cluster interior, leading to a gas 
inflow, at a rate M = Airr'^piVi. 



2.2 Stationary Structures, Stability and Observational Im- 
plications 

With the boundary conditions specified, Equations (4) and (5) can be inte- 
grated. Two example structures of clusters with different F2 are shown in 
Figure 1. For F2 = 3, the gas inertia is important in the cluster core; for 
F2 = 9, dark matter dominates throughout the cluster. In both cases, the 
temperature has a depression at about 10~^ of the cluster radius and in- 
creases at larger radii. For F2 = 9, the temperature first reaches a peak and 
then has a down-turn at the cluster out-skirt. The temperature down-turn 
at the cluster out-skirt is less obvious in the cluster with F2 = 3. Figure 2 
shows the X-ray surface brightnesses of the clusters. A bright hard X-ray 



core is present for F2 = 3, but is absent for F2 = 9. Figure 3 shows the local 
Jean lengths. For F2 = 2 both dark matter and gas are stable under the 
Jean's criterion. When F2 = 9, gas is stable throughout the cluster but dark 
matter can be unstable at the cluster core region. This implies that gas- 
free dark-matter collapse may occur under certain structure perturbations 
(Saxton and Wu 2008). 

To investigate the global stability properties of the cluster, we consider 
perturbations to the stationary structure and carry out a linear analysis. 
The Lagrangian perturbation is used, in which the hydrodynamic variables 
take the form: 



With this perturbative expansion, we obtain a set of differential equations, 
which are solved by direct integration after appropriate boundary conditions 
are imposed. The dimensionless eigenvalue of the dynamic equations is given 



The eigen-plane plots show sequential modes. The imaginary part of the 
eigenvalue {Si) alternates between positive and negative values. These pos- 
itive and negative 61 correspond to the growth and decay modes of the dis- 
turbances in the cluster. Our calculations show a trend that the strengths 
of the modes increase with the mode orders, implying that the higher-order 
modes grow or decay faster than the lower-order modes. 

The presence of radiative cooling of gas in galaxy clusters are expected 
to cause cooling gas inflow (Cowie and Binney 1977; Fabian, Nulscn and 
Canizares 1984). However, this phenomenon has not been confirmed by ob- 
servations. It has been speculated that the gas inflow may be halted as the 
radiative cooling may be compensated by certain efficient heating mecha- 
nisms. One argument is that hot electrons are injected into the intracluster 
medium by some means, e.g. jets from AGNs residing in the cluster (Mc- 
Namara et al. 2005). However, this would require some fine tuning. If the 
AGN injection is too strong the intra-cluster gas might be blown away; if 
the AGN injection is too weak it might be sufficient to quench the cooling 
flow. An alternative is the excitation of a cluster tsunami which could heat 
up the cool cluster core (Fujita, Suzuki and Wada 2004; Fujita et al. 2005). 
In the original tsunami model, the baryonic gas is trapped in a potential well 
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created by a static dark matter halo. The sound speed of gas is high in hot 
outer cluster region than is much lower in the cool cluster core. Disturbance 
(waves) generated in the outer cluster region may be subsonic initially, but 
when it propagates into the cooler cluster core, they could become super- 
sonic. The waves will therefore pile up at the cool core region, leading to 
wave tsunami and shocks. The dissipation in the shocks presumably causes 
heating and quench the cooling flow. 

Our analysis provides a more self-consistent cluster tsunami scenario, in 
that the dark matter halo provides a gravitational potential to trap the gas 
as well as acts a oscillator that drives the waves in the gas. The piling up 
of waves would occur as in the original tsunami model, but the tsunami is 
forced by the dark matter oscillations instead of self-cxcited. Our finding 
that higher-order modes grow faster implies that long wave-length oscilla- 
tions will break down into shorter wave-length oscillations. Thus, energy 
in the global cluster low frequency oscillation will be transfered to smaller 
scale oscillations, causing tsunamis and shocks which dissipate and heat the 
cluster core efficiently and evenly. 

3 Cosmic Walls 

3.1 Hydro dynamic Model 

Cosmic walls are semi-infinite two-component self-gravitating sheets. For 
self-interacting dark matter and baryonic gas, the equation of the state is 
Pi = Sip^' = piaf (here i = 1 for gas and i = 2 for dark matter). Since the 
dynamics are dominated by the dark matter, we may ignore the inertia of 
the baryonic component and simplify the conservation equations: 



— Vpi = -V(j); 



(8) 



Pi 



—P2 + V • P2V2 = ; 



(9) 



— V2 + (V2 • V)V2 Vp2 = -V0 . 

Ot P2 



(10) 



The Poisson equation is given by 




(11) 



where pi are the deviations from the equilibrium background density po- 



3.2 Soliton Formation 



Cosmic walls are sheet-like, and we use a quasi- ID treatment, such that 
d/dx » d/dy > d/dz = 0. We expand the density, velocity and potential 
in terms of a perturbative parameter e, i.e. 

p = pW+^^yj); (12) 
j 

V = ^£j7;(j) ; (13) 

j 

ct> = X]£V^^\ (14) 
j 

and consider stretched co-ordinates A = e^/^(x— Cgt), and r = e^^^t. Without 
loss of generality, we assume an isothermal gas component and set 71 = 1. 
The zeroth order expansion gives the effective sound speed 



Cs = ^72S2(p(°))^^-^+Ca2 , (15) 

where ^ = p^^"^ / {po — p^^^ ) . The first order non-linear expansion yields 
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Rescaling the variables leads to the KdV (Korteweg de-Vries) Equation: 

where the parameter p, is determined by the thermodynamic properties of 
the dark matter and the gas. 

The quasi-lD approximation can be relaxed by introduction of an addi- 
tional transverse co-ordinate y and its associated perturbation 77 = sy. The 
perturbation of the according velocity component would take the form: 

Vy = E^^^4^^- (19) 



With rescaling of the variables, the first order non-hnear expansion yields 
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which is the KP (Kodonistev-Petviashvili) Equation. The parameters H and 
p are determined by the thermodynamics of the dark matter and the gas. 

Note that the isothermal gas assumption, i.e. 71 = 1, is not essential for 
the derivation of the KdV and the KP Equations. The two equations can 
be obtained using a polytropic equation of state for the gas (Younsi 2008) . 

3.3 Astrophysical Implications 

A family of solutions to the KdV Equations and the KP Equations are the 
solitary waves (see e.g. Belashov and Vladimirov 2005). Solitary waves are 
non-linear waves, and they interact non-linearly. Linear superposition pre- 
dicts an amplitTidc doubling when two linear waves with equal amplitudes 
interact. However, collision of two solitons with equal amplitudes could 
form a metastable resonant pulse whose amplitude is four times that of the 
incident solitons. This non-linear superposition has important implications 
in the dynamics and structure formation of large-scale astrophysical sys- 
tems. For instance, when the resonant soliton resulting from the collisions 
of dark-matter solitions of positive densities causes a depression in the local 
gravitational potential, if the resonant time-scale is larger than the dynam- 
ical time of the baryonic gas, the gas will be trapped into this depression. 
Once the Jeans instability threshold is reached, the trapped gas will collapse 
and condense into visible structures. As there is a larger amplitude enhance- 
ment in the soliton interactions, it increases the efficiency of gas trapping 
and hence the subsequent condensation. 

In the quasi-lD case, we obtain a dispersion relation as 

It follows that the Jeans wave number 
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which is always positive for 71 > 0. Prom the dispersion relation, one 
can see that the system is stable for < A; < (47rGp^°)/^cri)^/^ and 



for k > kj, marginally stable for A: = and k = kj, but unstable for 
{4irGp^^^ /^ai)^/"^ < k < kj. As 71 >, bright solitons are formed when 
k > kj but dark solitons are formed for k < {iirGp^^^ /^af)^^'^ (Younsi 
2008). The bright solitons correspond to dark matter density pulses with a 
positive amplitude above the equilibrium background density, and the dark 
solitons correspond to dark matter density pulses with a negative amplitude. 
The two types of solitons can co-exist locally in the presence of high-order 
perturbations, leading to the formation of 'double layers' (Younsi 2008; see 
also Shukla and Mamun 2002). The formation of weak filaments (galaxies 
and galaxy clusters) deep inside a cosmic void has been a challenging is- 
sue in astrophysics. Conventional wisdom predicts no growth of structure 
inside large scale low-density regions, contrary to observation. The forma- 
tion of 'double layer' soliary waves in two-component self-graviating sheets 
may provide a viable model for the formation of such filament and void 
structures. 

4 Summary 

We construct a two-component model for large-scale self-gravitating sys- 
tems. The two components are self-interacting. Two example cases, galaxy 
clusters and cosmic walls, are presented. Our key findings are summarised 
as follows. When the dark-matter has a large internal degree of freedom, it 
can be Jeans unstable locally, implying the possibility of collapse under ex- 
ternal perturbations. Dark matter oscillations in galaxy clusters can lead to 
gas tsunamis in cluster cores, leading to heating of the intracluster medium. 
Dark matter solitons can be formed easily in cosmic walls, and resonance 
in soliton collision enhances gas trapping and condensation. Negative am- 
plitude solitons tend to be larger in size than positive amplitude solitons. 
The locally co-existence of the two types of solitons implies the possibility 
of filament formation in cosmic voids 
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DISCUSSION 

WOLFGANG KUNDT: Are you familiar with Tsallis Entropy? It is a 
generalisation of the entropy of the ordinary entropy of non-additive systems, 
which is called "non-extensive" by him. 

KINWAH WU: Yes, I know Tsallis' work on entropy and thermodynam- 
ics. Wc can incorporate part of Tsallis' idea through the assigned degree of 
freedom in the polytropic equation of state used in our formulation. 

JIM BEALL: Can you expand on your perturbative analysis of the system? 
Will the presence of an AGN jet change your calculations? 

KINWAH WU: It is possible to expand the perturbative analysis. What 
we need is to consider higher-order terms and keep the relevant equations. 

As for the presence of AGN jets, it will not change the formulation. We can 
treat the AGN jets as an energy source field. The energy injected by the jets 
will change the thermal structures of the galaxy clusters, which is certain. 




Figure 1: Stationary mass density and temperature profiles of model clusters 
witliF2 = 3 and 9 (left and right columns respectively). The mass inflow rate 
M = 10 MQyr"-'^ and the temperature at the outer boundary = 1 keV. Solid 
black lines correspond to dark-matter profiles and grey lines, gas profiles. The 
dotted lines indicate the cosmic mean density. 




Figure 2: X-ray surface brightness profiles of the two clusters in Fig. 1. Black lines 
indicate emission at the 0.1 — 2.4 keV band, and gray lines, the 2—10 keV band. 




Figure 3: The ratio of the radius to the local Jean length for the clusters in Fig. 1. 
Black lines indicate dark matter, and gray lines, gas. Arrows with labels 7, 8 and 
9 mark the marginally stable halos enclosing dark-matter mass of 10^, 10® and 
10^ Mr. 



Figure 4: Collision of two dark-matter solitons in the quasi-lD two-component 
cosmic wall. The two dark-matter solitons have equal amplitude (height) initially. 
They collide to form a resonant soliton which has an amplitude four times that of 
the incident solitons. This resonant dark-matter soliton has an infinite life-time. As 
the dynamical time of the baryonic gas, i.e. r (47rG'pi)~^/^, is finite, enhanced gas 
trapped and condensation will result from the formation of such resonant solitons. 



